arXiv:1501.07557v3 [cond-mat.str-el] 12 Jun 2015 


Excited States in Spin Chains from Conformal Blocks 
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We develop a method of constructing excited states in one dimensional spin chains which are de¬ 
rived from the SU{2)i Wess-Zumino-Witten Conformal Field Theory (CFT) using a parent Hamilto¬ 
nian approach. The resulting systems are equivalent to the Haldane-Shastry model. In our ansatz, 
correlation functions between primary fields correspond to the ground state of the spin system, 
whereas excited states are obtained by insertion of descendant fields. Our construction is based 
on the current algebra of the CFT and emphasizes the close relation between the spectrum of the 
spin system and the underlying CFT. This general structure might imply that the method could be 
applied to a wider range of model systems. 
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I. INTRODUCTION 

The description of quantum many body systems is an 
intrinsically hard problem due to the dimensionality of 
the Hilbert space, which depends exponentially on the 
system size. Despite the impressive success of numerical 
and approximate tech niqu es like Exact DiagonalizatiorP, 
Quantum Monte CarlcPSl^ and the Density Matrix Renor¬ 
malization GroupP, exactly solvable models remain es¬ 
sential in studying the physics of quantum many body 
systems. An analytical solution is not only indispens¬ 
able for benchmarking approximate methods, it can also 
elucidate the general structure of solutions and underly¬ 
ing physical principles. With the experimental advances 
in cooling and controlling the interactions of atoms, in¬ 
teresting model systems that have a closed theoretical 
solution may even be engineered and studied in the lab¬ 
oratory. 

Given the complexity of a generic quantum many body 
problem, the direct construction of exactly solvable mod¬ 
els appears as an appealing approach. In the past years, 
1-1-1 dimensional Gonformal Eield Theory (GET) has 
proven to be a powerful tool to construct model Hamilto¬ 
nians and corresponding ground state wave functions for 
continuum and lattice quantum Hall model^^^^. In this 
approach, a CFT is taken as the starting point to derive 
a quantum many body system. It can thus contribute 
to a systematic understanding of systems that admit a 
description in terms of a scale invariant theory, such as 
critical systems or systems with topological order that 
have a gapless edge spectrum. Theories for which such an 
analysis was carried out i nclude the SU{2)i^, S'0(n)il^, 
[/(l)qE^, and Wess-Zumino-Witten (WZW) 

theories. Another correspondence between CFT and a 
class of continuum and lattice quantum systems in two 
dimensions was studied in Ref. m In these models, the 
ground state exhibits a z = 2 Lifshitz scale invariance. 

In this paper, we describe a method of constructing ex¬ 
cited states of a spin system that is derived from a CFT. 
Correlation functions of fields in the CFT are interpreted 
as wave functions of states in the spin system. Following 


earlier studie^^l^^, we consider the SU{2)i WZW model 
and define a parent Hamiltonian for the state that corre¬ 
sponds to the correlator of N primary helds. In the case 
of periodic boundary conditions in one dimension, which 
we focus on here, the resultin g spin system is equivalent 
to the Haldane-Shastry modek^*^. 

Based on solutions to the Calogero-Sutherland 
modeP^lt^, excited states of the Haldane-Shastry Hamil¬ 
tonian were fi rst constructed as po lynom ials in the parti¬ 
cle basi j^^ l ^U It was realized lateJ^^I^ that these states 
are the highest weight states of the Yangian algebra, a 
hidden symmetry of the Haldane-Shastry Hamiltonian. 
Given the close relation between the Haldane-Shastry 
model and the SU{2)i WZW CFT, it was conjectured in 
Ref. [22] that there should be a correspondence between 
the states in the CFT and those of the spin chain. Here 
we show that this is indeed the case. Since our ansatz is 
based on the SU{2) currents, it is manifestly SU{2) in¬ 
variant. We relate our results both to the Yangian high- 
est weight states and to the spinon basis of CFT states^. 

Our ansatz mirrors the structure of the CFT, with the 
ground state corresponding to the CFT vacuum and the 
excited states corresponding to CFT descendant states. 
More precisely, we construct the excited states from the 
ground state by insertion of current operator modes into 
correlation functions of primary fields. We show that 
the Hamiltonian of the spin system is block diagonal in 
these states, with each block corresponding to a hxed 
number of inserted current operators. This allows us to 
obtain excited states by successively adding current oper¬ 
ators and diagonalizing the blocks. As the Hamiltonian, 
the Yangian is closed in the subspaces of states with a 
certain number of current operators. This implies that 
one can construct representations of the Yangian algebra 
in these subspaces. We explicitly construct the highest 
weight states for the analytically obtained eigenstates. 
Our ansatz for the excited states and its relation to de¬ 
scendant states of the CFT is summarized in Table |T| 

We perform the diagonalization of the Hamiltonian 
analytically for up to eight current operators, construct 
Yangian highest weight states in the obtained eigenstates, 
and confirm numerically for small system sizes that the 
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TABLE I. Structure of the states in the CFT and the spin sys¬ 
tem. The ground state of the spin system is constructed from 
the product of N primary fields 'l>s(z) = (j>si{zi) ... 
and corresponds to the CFT vacuum |0). Excited states are 
constructed from the ground state by insertion of current op¬ 
erator modes 7“^ ... J“\. 


CFT 


Spin system 


Ground 

state 


| 0 ) 


« (0|$s(z)|0) 


Excited 

states 




complete spectrum can be obtained in this way. We show 
that this construction can be done both for an even and 
for an odd number of spins in the chain. 

This paper is structured as follows: In Section |nj we 
review some properties of the SU{2)i WZW model and 
describe the states that form the basis of our construc¬ 
tion. We introduce the parent Hamiltonian in Section 


III and construct excited states, both analytically (Sec¬ 


tion IIIC) and numerically (Section IIIEl. We conclude 
in Section IIVI 


II. STATES FROM CONFORMAL BLOCKS 

In this section, we briefly review some properties of the 
SU{2)i WZW model and describe the correspondence be¬ 
tween conformal blocks and spin system wave functions. 

In addition to the identity, the model has one primary 
field with scaling dimension h = 1/4, the vertex operator 
4>s{z)- It can be constructed from the chiral part ^p{z) of 
a free, massless boson as 


algebrgP^ 

“ i^abcJm+n ( 3 ) 

with Eabc being the Levi-Civita symbol and 5ab the Kro- 
necker delta. 

A primary field (jjsiz) transforms covariantly with re¬ 
spect to conformal and SU(2) transformations. The 
latter is expressed by the operator product expansion 
(OPE) between the current J°‘{z) and (j)s{z'P^, 

r{z)(j)s{w )— XI 

s' 

Here are the SU{2) spin operators. They are related 
to the Pauli matrices cr“ by = (t“/2. 

The modes of the current operator give rise to a tower 
of descendant states 

(J“^..J“\)(0)|0), (5) 

with |0) being the CFT vacuum. States of this form build 
up the spectrum of the CFTI^. Note that it suffices to 
consider states obtained from the n = —1 mode of J°'{z). 
This is so because one can successively rewrite a higher 
order mode n > 0, in terms of the lower order modes 
J-n+i by means of the Kac-Moody algebra (cf. 

Eq. (§), 

= ( 6 ) 

In this work, we show that the spectrum of the spin 
systems organizes in the same way in terms of current 
operators. The excited states we obtain are linear combi¬ 
nations of conformal blocks containing current operator 
modes JO*) ... J“\. We give a summary of these states 
in Table U and describe the construction of states from 
conformal blocks in the next subsections. 


Here s = ±1 corresponds to the two components of the 
vertex operator, and the colons denote normal ordering. 
(The holomorphic field />s(z) has an anti-holomorphic 
counterpart (paiz). In the following we only consider the 
holomorphic sector.) The value q G {0,1} corresponds to 
the two sectors of the CFT: g = 0 if the operator </s(z) 
acts on a state that has an even number of /i = 1/4 pri¬ 
mary fields and q = 1 for a state with an odd number, 
respectively. 

In addition to conformal invariance, the SU{2)i WZW 
model has an SU{2) symmetry, which is generated by the 
current operator J°‘{z). Its Laurent expansion defines 
modes J“, 


OO 

J“(z)= ^ (2) 

n— — oo 


where a G {x,y,z}. They satisfy the Kac-Moody 


TABLE 11. Summary of the different towers of states ob¬ 
tained by insertion of current operator modes, 'l's(z) denotes 
the product of N primary fields, 'I>s(z) = cjAsi (-^i) • ■ • 4>siv{zn)- 
Using the OPE between (j)s{z) and J‘^(z), the wave functions 
for these states can be written as the application of k Fourier 
transformed spin operators to the state without current op¬ 
erators, as explained in Section EE 

Tower of states See Eq. 


A odd {4-,(z)(J“^.. J“\0,J(O)) @ 


(14 


N even (4-s(z)(J!"i ... J“\)(0)) _ 

A even {cl)s^{oo)^siz){Jl\ ■ ■ ■ (^ 

A odd (0,„„(oo)4-s(z)(J/"i... J“\)(0)) ^ 


A. State Obtained from a String of Vertex 
Operators 

We consider an even number of vertex operators (/g. (zi) 
{i = 1,...,A) that each transform under a represen- 
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tation of SU{2) generated by spin operators tf, as ex¬ 
pressed by the OPE of Eq. Q . The key idea is to view 
the correlation function of vertex operators in the GET 
as the wave function of a system of spin 1/2 degrees of 
freedom on a lattice, 

IV'o)= ^ iI}o{si,...,sn)\si,...,sn), (7) 

SI...SN 

where , sn) is given by 

V'o(si,---,SAr) = ((/si(zi)...(/s„(zAr)). (8) 

Here, 

(j)^^ (Zj) = (g) 

and |si,...,SAr) with Si = ±1 is the tensor product of 
eigenstates |si) of the z-component of the spin operator 
tf. The coordinates Zi define the lattice positions in the 
complex plane and are kept fixed. In contrast to the 
continuum case, there is no spatial degree of freedom 
in the basis states |si,..., s^v) and therefore no integral 
over the positions in Eq. Q. This is why we use the 
notation '0o(si) ■ • ■ sn) without the coordinates Zi for the 
spin wave function. 

The correlation function of N vertex operators is given 
bjl^ 

V'o(si, ■ ■ ■ ,SAr) = (10) 

N 

i<j 


B. States Obtained from Vertex and Current 
Operators 

In analogy to the tower of GET descendant states 
(J“\ ... J“\)(0)|0), we define a tower of spin states by 
insertion of modes of the current operator J“(z) into a 
correlation function of vertex operators, 

|V'afc...ai) = E (si,..., sjv)|si,..., Sn), 

SI...SN 

where (si,..., sn) is given by 

...oi (si , . . . , Sn ) 

= (</.,(zi). ..<PsA^N){r\ ■ ■ ■ r\)(0)). (14) 

It is possible to obtain the states ipak...ai by applying 
spin operators to tpo- To see this, we first note that 

(4>.(z)(J^i?)(0)) = —(4>.(z)J“HH(0)), (15) 

27rz Jo w 

with $s(z) = (/s^(zi).. .4)sj^{zn) and B being an arbi¬ 
trary operator. Applying this relation to the definition 
of Tpak-.-ai we obtain 

Yak...ai (si, • . ■ , Sn) 

= (zi).../),„(z^)(J“^..J“\)(0)) 

= i^<f —{M^)rHw){rr...r\)m 

27rz Jq w 

T ^ r A 

j=i 

(16) 


where Ss is 1 if = 0 and 0 otherwise. Xs is the 

Marshall sign factor, 

N 

( 11 ) 

p^l 

which ensures that the state 'tpo is a spin singletf^El^ 

AT 

= ( 12 ) 

2=1 


Using the OPE between a current operator and a primary 
field (cf. Eq. @)> we get 



Successive application of the same argument results in 


Note that the condition requires the number of pri¬ 
mary fields (pSjYj) llio correlator to be even. 

The positions Zj can, in principle, assume any value on 
the complex plane. In the following we mostly consider 
N spins uniformly distributed on the circle, 

Zj=z\ withz = e2’^*/^. (13) 

In this case, the state ipo has a momentum of tt if N/2 is 
odd and 0 if A^/2 is everf^, see also Appendix [d| 




afc...ai 




(18) 


If the positions Zj are uniformly distributed on the circle, 
we can express this result in terms of Fourier transformed 
spin operators u/. 


u 


a 

I 


N 


f^a^27TijllN 


(19) 
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With z = in Eq. (18), we have 


V'afc...ai = .. .u“\'0o- (20) 

Therefore, each additional insertion of J°Ly changes the 
momentum of the state by 2ttIN. 


C. States with Additional Spins at Zero and 
Infinity 

We define an additional class of states by insert¬ 
ing two extra vertex operators into the correlator 
{4>si ( 2 : 1 )... {^n)), one at 2 = 0 and one at 2 = 00 , 

l^so.soo) ^ ^ (21) 

Si,...,SAr 

V'o“’'*”(si,...,SAr) 

= {oo)(l)s, ( 21 ) . . . {zn) 4 ‘so (0))- 


This ansatz thus corresponds to the tower of CFT de¬ 
scendant states (J“\ ... J“\(/)s)(0)|0). 

An argument similar to the one given for ipat-.-ai shows 
that 




SO.Soo 
ak -ai 



N 


+ E 



X 





So,S, 

0 


(24) 


Note that the terms /Zoo do not contribute in the limit 
Zoo —>■ 00 . On the unit circle, we have 

(25) 

in terms of Fourier transformed spin operators. 


In the Riemann sphere picture, the additional spins are 
added at the south and north pole, respectively, while 
the N spins at the unit circle are located at the equator 
(see Fig. [^. 


N 

I 



D. Odd Number of Spins 

A correlation function of vertex operators 
{(psiizi) ■ ■. ips^iz^N)) is only non-zero if the sum 
'S* vanishes (cf. Eq. ©)■ This property, the 
charge neutrality condition, implies that the vertex 
operators in the correlator need to have a net charge of 
zero. As a consequence, the number of vertex operators, 
and therefore the number of spins, needs to be even. 
We can, however, still consider a model with an odd 
number of spins at the circle, by adding an extra vertex 
operator that compensates the excess charge at the 
circle. Inserting an additional vertex operator at 2 = 00 , 
we obtain the state 


FIG. 1. (Color online) The Riemann sphere with N spins at 
the equator (unit circle) and two additional spins, one at the 
north pole {z = 0) and one at the south pole {z = 00 ). 

The wave function is 

N 

(si,..., sat) cx n 

n—1 
N 

X Y[iZn-Zmr-^-^\ ( 22 ) 

n<.m 

where Jg = 1 for sq + Soo + Sj = 0 and dg = 0 
otherwise. 

Note that the extra fields inserted at zero and infinity 
are primary fields. If we insert additional current oper¬ 
ators, we generate descendant states. We thus define a 
tower of states on top of j 

V'a“’.^.ai(si,---,S7V) (23) 

= (</>S3o(oo)<l^si(^l) • ■ •<(>Siv(^Ar)('^-l ■ • • -^-\<(>so)(0))- 


■!/'o”(si,...,SAr) 

= {oo)(l)si ( 21 ) . . . (zn)) 

N 

(X 5g(-l)(*“+i)/2^g Hiz, - (26) 

i<i 

where (5g = 1 if s^o + Sti = 0 and 5g = 0 otherwise. 
The wave function has spin 1/2, 

3 iv 

, with T“ = ^ (27) 

i=l 

This is a consequence of V'o” being a singlet of the total 
spin including the point 2 = oo, (t^ -I- T“)'!/g“ = 0. As 
in the case of an even number of spins, we define a tower 
of states by insertion of current operators, 

V’ar..ai(si.---,SAr) 

= (<(>s„„(oo)(/si(2i) . . . />g„(2Ar)(J“'i . . . T“\)(0)). 

(28) 
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We obtain a second class of states by inserting the addi¬ 
tional vertex operator at z = 0 instead of z = oo, 

N N 

oc ^sXs n W{z, - (29) 

i=l i<j 

and the corresponding tower 

= (</)«, (zi). (z^)(J“\ ... r\K)m- (30) 

On the unit circle, the two towers of states can be 
written as 


and 




dj^ ,... ,<21 


(si,..., sn) = • ■ • U-\ipo 


(31) 




<2fc,...,a 


^(si,...,SAr) =u°L\---u°L\Yi 


(32) 


with 


Wij = 


Zj -I- Zj 
Zi - Zj ' 


(35) 


We have used the notation J2i{^j) ^ sum over all 

ie {l,...,7V}\{j}. 

This allows for the definition of a parent Hamiltonian 
H of 




(36) 


Note that H is positive semidefinite and ipo is an eigen¬ 
state of H with zero energy. 

It is knowiP that H is closely related to the Haldane- 
Shastry Hamiltonian Hhs j if the spins are uniformly dis¬ 
tributed on the circle, 


Hks = 


ja ja 

‘'i 0 


9 ^ . 2 

sm 


/ (i-j)TT \ 
\ ^ ) 


= H 


N + 1 
6 


T‘^T‘^ + Eo. 


(37) 


respectively. 

Comparing the wave functions of Eq. 


(26) and Eq. 


(29) with the case of an even number of spins (Eq. (10) 
and (22)), we conclude that Yo 
ipZ of ■ 00 °*“ j respectively. 


is the analog of Yo and 


III. PARENT HAMILTONIAN AND 
SPECTRUM 


Here T“ = is the total spin and Eq = —{N^ + 

57V)/24 the ground state energy of the Haldane-Shastry 
Hamiltonian. Note that tjjo is annihilated by E[ and also 
by T“T“ in Eq. (37) since it is a singlet. Therefore, ipo 


is the ground state of the Haldane-Shastry Hamiltonian. 
From now on, we will work with the Hamiltonian 


H = H- 


N +l 
6 


rj^CLrj^a. 


(38) 


In this section, we introduce the parent Hamiltonian of 
ipo, discuss its equivalence to the Haldane-Shastry model 
and construct its spectrum from GET current operators. 
Furthermore, we relate our ansatz to the multiplets of the 
Yangian algebra and the spinon construction of the GET 
Hilbert space. From now on we assume a uniform, one¬ 
dimensional lattice with periodic boundary conditions, 


Zj = 


(33) 


with j G {1,...,7V}. 

A. Parent Hamiltonian 

In earlier work on the Haldane-Shastry modeP^, an 
operator C“ was constructed that annihilates the wave 
function ipo^ CfYo = 0- This operator was later obtained 
in a more general setting from the SU{2)i WZW model 
using null vector^. In terms of spin operators, 

^ E + i£abct\tf) , (34) 


dropping the constant Eq. 


B. Block Diagonal Form of the Hamiltonian 


We now systematically construct excited states of % 
from conformal correlation functions. Specifically, we 
build the excited states as linear combinations of the 
states '0afc...ai (cf. Eq. (10) and ([I^). 

The key to this construction is that the Hamiltonian 
does not couple states with a fixed number of current 
operators to states with a different number of current 
operators, i.e. the Hamiltonian is block-diagonal in this 
basis. 

Therefore, we can diagonalize the Hamiltonian in the 
subspaces of states with a certain number of current op¬ 
erators. It is not necessary to construct the Hamiltonian 
in the full Hilbert space of dimension 2^ in order to find 
eigenstates beyond the ground state. Rather, we obtain 
eigenstates by successively adding current operators and 
diagonalizing the blocks. 


Let us now show that Hipa 


is a linear combina¬ 


tion of states obtained from ipo by insertion of k current 
operators. 
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Recall that the insertion of k current operator modes 
J°l\ {j = 1,..., fc) into the correlation function of ver¬ 
tex operators is equivalent to the successive application 
of Fourier transformed spin operators u°^i to the ground 
state V'o (cf- Eq- (20l). Therefore, we have computed 
the commutator between % and u°ii by an explicit ex¬ 
pansion of the Hamiltonian in terms of Fourier modes 
(cf. Section 1^ of the Appendix). The result of this 
calculation is 


N 


3C? 


[H, i] ={N- l)u“ 1 ^ ^ + ieabcu\T-. (39) 

Z Zi 


Combining Eq. (39) and Eq. (41), we obtain 


T~L'4’ak...ai 

k 

= ^ [n, uZi] ^“7^ 


r—1 


q-1 


2<q<r<kn—0 

+ E ■ ..a-r+iaqar—i ...aq^iaraq — i...ai 

l<q<r<k 

‘^^arttq V'a,.. .(lrJ^±CClj- — l .. .Clq_|. 1 CCLq — \ . . .CL\ 


From this we can already conclude that the energy of a 
state with one current operator is iV — 1, 


-I- '0 

- '0 


dfj,.. .d'P^XdqCL'pCL-p — x...dq^±dq — x...dl 


dk---dr+ldrdqdr_l...dq^ldq_i ...di 


HU010O = [n, u0i] ^o = {N- l)u0iV’o, (40) 


with 


-I- 2SN2Sar.a^'ip 


dk---d-r^ldr-l ...dqj^idq—i...di 



since H, Cf, and T° annihilate the ground state ipQ. 

We need to know how Cf acts on ipa^-.-ai to determine 
the energy of states with more than one current operator. 
As we show in Section of the Appendix, 


Ci fpak...ai 


q=l 


Zi 


d^~...dqj^ldq — i...di 




n+1 

q—2 n—0 

Tflfc J^q + l TC J^q-l 


X ($,(z)(j“^i... rrjfr\-^ ... j“\)(o)), (41) 


with 


T^qr,n 

dk...di 

= 2 {^,{z){r\. .. 

jf'-rv ...jf\m) 

- . .. Jir 

jf^rv ■ ■ ■ J-\m) 

+ 27V5„+2^ea.a,c($s(z)(J“^l . . . • • ■ J-t" 

( 45 ) 

In the last term, Sn +2 = 1 if (n -I- 2) mod A^ = 0 and 
6 n +2 = 0 otherwise. 

Let us now argue that this expression contains k cur¬ 
rent operators of order —1. There are two terms that are 
not yet explicitly written in the desired form, namely the 
term proportional to 5 n 2 and the term abbreviated by 

Tp<i'c,n 
■^dk...di • 

The operators (and respectively), can be 

written as a linear combination of n-|-2 current operators 
of order —1 by repeated application of Eq. (|^. On the 
other hand, the operators J“’' (and Jf, respectively), can 
be commuted to the right, using the current algebra (cf. 
Eq. (§) 


.. .(ps^izN), and 
(Kf), = ^{Sab-iSabctd- 


(42) 

(43) 


[J“, Jti] = iSabcJf-l + ^SabSn- 1 , 0 - (46) 

The resulting terms either have n current operators less 
or vanish because Jfn\ 0 ) = 0 for to > 0. 

The total number of current operators in the first two 
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terms of is therefore 

k — 2 + n + 2 —n = k. (47) 

J-l-2 

If fc > -/V, there can be a contribution from the third term 
in namely if n+2 = mN for m £ {1, 2,... }. This 

term can be written in terms of k—mN current operators. 
Note, however, that the space of states with k current op¬ 
erators contains the space of states with k — mN current 
operators. The reason for this is that a current operator 
Jti corresponds to the application of a Fourier trans¬ 
formed spin operator u°_i, for which u°ii = 

This argument also applies to the term that is propor¬ 
tional to 5 i^2 - 

Thus, the Hamiltonian is block diagonal in the states 
with a fixed number of current operators. 

Note that this observation does not follow from trans¬ 
lational invariance only. Translational invariance implies 
that the Hamiltonian does not mix states with different 
lattice momenta. Since each current operator in '0afc...ai 
contributes a unit of 2tt/N to the momentum, it follows 
from translational invariance that 'H'ipak...ai is a linear 
combination of states with k mod N current operators. 
The above considerations moreover show that it is pos¬ 
sible to write 'H'4’ak...ai as a linear combination of states 
with strictly k modes Jli- In particular, it is not nec¬ 
essary, to include terms with a higher number of current 
operators. This allows us to block-diagonalize the Hamil¬ 
tonian starting with the smaller blocks, i.e. those with a 
small number of current operators. 


C. Eigenstates from Current Operators 


We have solved the eigenvalue equation of Eq. (441 for 


up to eight current operator modes analytically. Since 
the Hamiltonian is SU(2) invariant, we have decomposed 
the eigenstates into different spin sectors. The momenta 
of the states are directly related to the number of cur¬ 
rent operators, with each current operator changing the 
momentum by 2'k/N. We summarize our results for up 
to four current operators in Table [lll| 

At level one we find a triplet (spin one), at level two 
a singlet and a triplet, at level three we find a singlet 
and two triplets with different energies. A spin two state 
appears at level 4 as the symmetric traceless part of a 
state with two non-contracted SU (2)-indices. 

Note that the number of eigenstates is smaller than 
the number of possible combinations we can build with 
k current operators. The reason is that some CFT states 
... J“\|0) are null, such that the norm of the corre¬ 
sponding spin state vanishes. At level three, for example, 
there is the null state - 2^’^“). 

The number of states that we find with a certain num¬ 
ber of current operators is in agreement with the charac¬ 
ters of the SU{2)i algebrE^ At level 0, 1, 2, 3, and 4 we 
find 1, 3, 4, 7, and 13 states, respectively. The size of the 


matrices that need to be diagonalized at a given level in 
current operators thus corresponds to the characters of 
SU{2),. 

When considering the action of the Hamiltonian on 
states with k current operators, there are two types of 
terms that depend on the number of spins N, cf Eq. 
(44). The first one, k{N — l)V^ofc...ai, is already diago¬ 


nal. All other terms only appear if fc > A^. These terms 
are strictly upper triangular in the sense that they can be 
written in terms of k—mN current operators with m > 0. 
This upper triangular structure is preserved by a diago- 
nalization of all other terms. Therefore, only the diagonal 
term k{N — l)V'afc...oi contributes to the A^-dependence 
of the energies. 

We thus arrive at an A^-independent representation of 
the energies by subtracting the contribution k{N — 1) 
from the energies, 


.> 0 , 

10, k = Q. 


(48) 


(We have also rescaled the energies by 1/k for conve¬ 
nience.) 

The shifted and rescaled energies E as well as the spin 
content of the corresponding eigenspaces are plotted for 
up to eight current operators in Fig. 

Note that for a given N, some of these energies do not 
occur because the corresponding states have zero norm. 
These null states appear dependent on N and in addition 
to null states identified at the CFT level. The occurrence 
of additional null states reflects the fact that the CFT 
Hilbert space is infinite, while the spin system’s Hilbert 
space is hnite. 

A necessary condition for a particular state to be non¬ 
null is given by i? > 0. The region of energies, where 
A > 0 is satisfied is indicated by the bands in Fig. 
Depending on N, we find certain states that violate this 
condition and are therefore null. These states then lead 
to further null states at higher levels through the appli¬ 
cation of additional current operators. We hnd that if a 
state corresponding to a certain energy is null for a num¬ 
ber of spins N', it is also null for N spins with N < N'. 
For each energy level, the highest N' that we found ex¬ 
ploiting the energy condition F > 0 is given in Fig. 
below the spin content. 

Let us give an example for the notation used in Fig. 

At k = 7 and E = —6 we find two spin 1 states and 
one spin 2 state. One of the spin 1 states is null for all 
N < 6, the other spin 1 state and the spin 2 state are 
null for all N <8. This multiplet is shown in Fig. [^as 


-6 

1 © 1 0 2 

6 8 8 


(49) 


Note that the violation of the energy condition E > 0 
is sufficient for a state to be null. A complete separation 






TABLE III. Eigenstates of H in terms of states obtained by insertion of k current operator modes for A: < 4. The momentum 
of the ground state is po = tt if N/2 is odd and po = 0 if N/2 is even. 


k State 

Null for Energy 

Spin Momentum Number of states 

o 

o 

II 

o 


0 

0 

Po 

1 

1 = Ipa 


Y- 1 

1 

po - 2'k/N 

3 

2 “ 3(5jv2'0o 

N <2 

2(Y- 3) 

0 

Pa — 47r/Y 

1 

2 - 'y ^(^cd'^cd 

N <2 

2(Y- 3) 

1 

Po - 47r/Y 

3 

3 ^ ^cde'^cde 

iV < 4 

3(Y- 5) 

0 

Po — Gn/N 

1 

3 ^ ^'4^cac “h V^cca) 4(5iV2'0a 

Y < 4 

3(Y- 5) 

1 

Po — Gn/N 

3 

3 (pi®' = - V’acc) + 25^200 

N <2 

3(Y- 3) 

1 

Po — Gn/N 

3 

4 = I]cd(50cdcd - 30cddc - 20ccdd) -b 165iV40o -b 12(5iV20o < 6 

4(Y - 7) 

0 

Po - 8tv/N 

1 

4 <^a —■ (4£acd’0cdee ^^acd'^cede.^ 

N <6 

4(Y - 7) 

1 

Po — Sn/N 

3 

4 (pi®) = "^^^{ipcddc — ’tpccdd) “ 125jV4'0O -b 6^iV2'0O 

Y < 4 

4Y - 18 

0 

Po - 8tv/N 

1 

4 ■ y vgjjg (^acd'0cdee 4“ ^^acd'^cede') 

Y < 4 

4Y- 18 

1 

Po - 87V/N 

3 

4 Xjc ( 2 ~h V^bacc) '^^ab '4^ddcc'j 

N <2 

4Y- 10 

2 

Po — 871/N 

5 
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FIG. 2. (Color online) Analytically calculated energy levels and their spin content as a function of the number of current 
operators k. E is defined sls E = {E — (N — l)k)/k for k > 0 and E = E = G for fe = 0. The three rows in the boxes next 
to the energy level correspond to the calculated values of (1) the energy E (2) the spin content (3) a value for the number 
of spins N' so that the corresponding state is null for all systems with N < N' spins. The inequalities indicate for which N 
the condition E > 0 is satisfied. For a given number of spins N, all eigenvalues that lie in bands where the inequality is not 
satisfied correspond to null states. By applying additional current operators to these null states, we could identify further null 
states that do not violate the energy condition. 


of null states requires the computation of inner products 
between our ansatz states at the level of the spin system. 
As shown in Ref. [51 the spin correlation functions in '00 
can be computed by solving linear algebraic equations. 
Using these correlation functions, one could compute the 
norms of the states 0afc...ai numerically, even for large 
system sizes. 


D. Construction of Yangian Highest Weight States 
and Comparison to Spinon Basis 


The spectrum of the Haldane-Shastry model has previ¬ 
ously been constructed by exploiting a hidden symmetry 
of the Hamiltonian, which is generated by the rapidity 
operator A"!^, 


A“ = ^ ^ w^j£abctitj = mjiu X (50) 
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Using = 0, it follows that 

o N 

A“= 4 E^^ ( 51 ) 

1=1 

The rapidity A“ and the total spin r“ form a basis of the 
Yangian algebra. They both commute with the Haldane- 
Shastry Hamiltonian, but the rapidity A“ does not com¬ 
mute with the total spin operator This is the rea¬ 

son for the degenerate energy levels formed by multiplets 
with different total spirP^. 

The key to the construction of eigenstates in this ap¬ 
proach is the notion of a Yangian highest weight state 
h, which is annihilated by A+ = A® -f iA^ and = 
-I- iTy. A multiplet of states with the same energy is 
then given by application of powers of A“ = A“ — iA^ to 
h. 

In order to relate our method to this approach, we have 
computed the action of A“ on ipak...ai using the decou¬ 
pling equation derived in section of the Appendix. We 
find 


2A Ipak-.-ai — 'y l)*£aa,c'*/’afe...a,+ica,_i...ai 

9=1 

k 

+ E ’^^aqO.c'^cak 
9=1 

+ 1 : 1 : or..... (52) 

q—2 n—0 

with 

^afc...aq ~ 

■ ■ ■ J-iKo)) 

- • ■ • ^-\)(o)) 

+ 2NSn+iieaa,Mz)iJl\ ■ ■ ■ 

(53) 


TABLE IV. Highest weight states of the Yangian algebra in 
terms of the eigenstates of the Hamiltonian given in Table |IH| 

k State Energy Spin 


0 



0 


0 

1 


+ 

N - 

1 

1 

2 


, • (3) 

2{N 

-3) 

1 

3 

(5) 

(pi 

1 • (5) 

3{N 

-5) 

1 

3 

(6) 

(pi 

, • (6) 

3{N 

-3) 

1 

4 

(8) 

(pi: 

, • (8) 

4(V 

-7) 

1 

4 

(10) , . (10) 

(pi: + *(Pa 

AN- 

- 18 

1 

4 

(11) , o (11) 1 O' (11) 
H” 2(^xa: H” 

AN- 

- 10 

2 


Note that our ansatz is manifestly SU{2) invariant, 
whereas the construction in terms of highest weight states 
is not. The states tpak...ai have a simple form in the sense 
that they are created from the ground state by successive 
application of Fourier transformed spin operators (cf. Eq. 

@). 

We have compared the highest weight states of the 
spin system to the highest weight states of the Yangian 
algebra in the CFT, which is spanned bjEH and 


m—1 


(55) 


We find that the highest weight states that we computed 
in the spin system are precisely the highest weight states 
of Jq and Q“, when seen as states of the CFT under the 
correspondence of Table |Tj 

Finally, let us relate our ansatz to the spinon basis 
of the CFT Hilbert spac^^. In this approach, states 
are not constructed by application of modes of the affine 
current </“„ to the CFT vacuum, but by application of 
modes 4>s,-m of the primary held 4>s{z). These modes 
are dehned with respect to the Laurent expansiorP^ 

<t)s{z) (56) 


Furthermore, we have (cf. Eq. (C61 in Appendix 


N 


T l/'afc...ai ) ^aa„c'0afc...a„+icao_i...gi ■ (54) 

9=1 


This shows that the Yangian, like the Hamiltonian, 
leaves the subspaces of states with a hxed number of cur¬ 
rent operators invariant. It is thus possible to write the 
highest weight states of the Yangian algebra in terms of 
the states ipak-.-ai- We have computed the highest weight 
states for up to four current operators and expanded the 
result in the eigenstates listed in Table [IIl| These states 
thus correspond to the eigenstates constructed by Hal¬ 
dane in Ref. nn which have a polynomial form in the 
particle basis and were identified as the highest weight 
states of the Yangian algebra in Ref. |2H We summarize 
our results in Table HVl 


Here q € {0,1} corresponds to the sector of the CET 
Hilbert space on which (j>s{z) is acting: q — 0 for a state 
built from an even number of spinon modes and < 7=1 
for a state with an odd number of spinon modes. The 
authors of Ref. [24] derived generalized commutation re¬ 
lations for the modes of J°‘{z) and (j>s{z), which can be 
used to rewrite states built from current operator modes 
in terms of spinon modes At level fc = 1, for ex¬ 
ample, the state ... (j)siv{zN)J-i{0)) corresponds 

to 

Xl(^“)s«'('('si(^l) (^iv) (0)). 

ss' 

(57) 

A general state with k current operators of order — 1 
will be a linear combination of terms with spinon modes 
^ai ,—mi ■ ■ ■ mi with k — 1 
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E. Complete Spectrum 


In the previous subsections, we have analytically diag¬ 


onalized the Hamiltonian in the states 'ipa 


for fc < 8 


and N even. For a given k, the size of the matrices that 
need to be diagonalized and therefore the complexity of 
the analytical calculation does not depend on the num¬ 
ber of spins N. However, it becomes increasingly diffi¬ 
cult for larger k. In order to test if our method yields 
all excited states or just a subset thereof, we have thus 
constructed the states ipak...ai numerically for small N 
and performed a numerical diagonalization of the Hamil¬ 
tonian in that basis. Our numerical calculations confirm 
that the complete spectrum is indeed obtained from our 
ansatz states. 

In an analogous way to we studied the tower of 
states obtained fro m ^^° ’^°° by insertion of current op¬ 
erators (cf. Section HC). We find numerically that the 
complete Hilbert space is generated, starting from 
with So,Soo G and successively inserting current 

operator modes. As for 'i/jg, the Hamiltonian is block 
diagonal in the states with a fixed number of current op¬ 
erators. 

Furthermore, we constructed the spectrum nume ricall y 
for the case of an odd number of spins (cf. Section HD). 
We find that the Hamiltonian is block diagonal in the 
states V’a“...ai and and that the complete Hilbert 

space can be constructed from states of this form. 

The numerically obtained spectra are shown for N = 7 
and N = 8 spins in Fig. 




k 


40 
gg30 

S20 

<D 

10 
0 
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FIG. 3. Numerically calculated spectra for the odd-A spin 
chain {N = 7, left panels) and the even-A spin chain (A = 8, 
right panels). The upper panels show data obtained by apply¬ 
ing current operators to the states and tpo, respectively. 
In the lower panels, the current operators were applied to 
the states and , respectively. The horizontal axis 

shows the number of current operators k that are present in 
the corresponding states. The numbers above the levels indi¬ 
cate their degeneracy. Notice that the spectrum is the same 
in the upper and lower plots, but a given state does not nec¬ 
essarily appear at the same number of current operators. 


ber of current operators have a larger energy. This means 
that by inserting few current operators, we get access to 
low-lying excited states. Note, however, that this relation 
does not hold in a strict sense; it may happen that, by 
inserting additional current operators, we obtain states 
with a lower energy. 

We observe that the spectra shown in the upper and 
in the lower panels of Fig. are the same. This means 
that for both N even and N odd, the complete spectrum 
is obtained for either of the two classes of states. 

The top right panel of Fig. corresponds to our ana¬ 
lytical results shown in Fig. Note that not all states 
of Fig. [^appear in the numerically calculated spectrum, 
because certain states of Fig. [^are null for N = 8. This 
is the case for the all states with E < 0, which appear in 
the bands below E = —7 in Fig. The remaining null 
states have the values 

k E E Spin 

6 -7 0 00 1 

7 -f 1 00 1 

7-6 7 102 

8 -6 8 00 1^0 2 

and can be obtained from null states violating the energy 
condition by the insertion of additional current operators. 

We have numerically computed the maximal number 
of current operator insertions of order —1 needed to gen¬ 
erate the complete spectrum as a function of A. The 
results are shown for the tower of states built on '0o (A 
even) in Fig. In this case, our numerical data suggests 
that (A/2)^ current operators are sufficient to obtain the 
complete spectrum. 

20 [- ^^^^^^- 

15 • 


5- 

gl-^^^^^^^- 

012345G789 

N 


FIG. 4. Maximal number of current operator insertions femax 
of order —1 needed to obtain the spectrum of the Hamiltonian 
dependent on the number of spins A. The data is consistent 
with kmax = (A/2)^. 


Finally, let us comment on how the first excited states 
of the Haldane-Shastry model for A even are related to 
the states 'tpak...ai- Except for A = 2, these cannot be 
given by the states at level one, ipai, which have E = 
A—1. The reason is that the states are eigenstates 

with a lower energy of A = A/2. This follows from 


We observe a tendency that states with a higher num- 




sq,s, 

0 


iKSUtl-4W 


(58) 
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and 

Note that the latter equality is only valid on the circle 
with Zj = Our numerical calculations indicate 

that the states are indeed the first excited states. 

Furthermore, we find a multiplet with spin content 0 © 1 
and energy N/2 at k = N/2 among the states constructed 
from ipak.-.ax-, both in the numerical and in the analytical 
calculations. These multiplets appear at E = —2k + 2 
in the analytically computed spectrum shown in Fig. 
Noting that the states can be decomposed into a 

singlet and a triplet, this suggests that the first excited 
states are given by linear combinations of states tpak-.-ai 
with k = N/2 current operator modes. 


IV. CONCLUSION 

We studied a model of a quantum spin chain that is 
constructed from the SU{2)i WZW model and is equiva¬ 
lent to the Haldane-Shastry model. We have shown that 
it is possible to construct excited states of the spin sys¬ 
tem from CFT current operators in a way, which directly 
reflects the structure of the CFT spectrum. 

In the approach pursued in this work, correlation func¬ 
tions of operators in the CFT are interpreted as wave 
functions of spin states. Based on earlier work, where 
the correlation function of primary operators were inter¬ 
preted as the ground state wave function, we provided 
a method of constructing excited states by inserting de¬ 
scendant fields. 

In the case of an even number of spins N, we have diag¬ 
onalized the Hamiltonian analytically for k < 8 inserted 
current operators. Depending on N, we identified cer¬ 
tain states that are null but correspond to non-null CFT 
states. These additional null states occur due to the fi¬ 
nite size of the spin system’s Hilbert space. Our method 
of detecting them is based on a sufficient criterion for a 
state to be null. In order to test if a given state is not 
null, one could use the known algebraic equations for spin 
correlation function^ in 't/o to compute the needed inner 
products numerically, even for large system sizes. 

We have given numerical evidence that this method 
yields the complete spectrum for a given number of spins 
N, independent of which primary state we build the 
tower states from. Furthermore, we have shown numeri¬ 
cally that a similar construction can be made for an odd 
number of spins. 

Our manifestly SU{2) invariant ansatz is compati¬ 
ble with the construction of eigenstates of the Haldane- 
Shastry model as multiplets of the Yangian algebra in 
the sense that the Yangian operator does not change the 
number of current operators. This allowed us to explic¬ 
itly relate the excited states constructed from current 
operators to the highest weight states of the Yangian op¬ 
erator. Furthermore, we have argued that our ansatz 


wave functions with k current operator modes of order 
— 1 can be rewritten in terms of a wave function with 
spinon modes whose mode numbers sum up to —k. 

In the case of the SU{2)i WZW model, which we stud¬ 
ied here, the resulting spin system is equivalent to the 
Haldane-Shastry model. Thus our method provides an 
alternative way of constructing the excited states of the 
Haldane-Shastry model, which emphasizes its close rela¬ 
tion to the underlying CFT. We expect that this method 
could be generalized to the SU{n)i WZW model, which 
is related to the SU{n) Haldane-Shastry spin chaiiP^. 
Another generalization of our construction could be pos¬ 
sible within the Laughlin lattice models with filling frac¬ 
tion l/cF. Even though these systems do not have a 
Yangian symmetry for q > 2, part of the spectrum is de¬ 
scribed by integer eigenvalues. We have carried out ex¬ 
emplary numerical calculations for g = 3 which indicate 
that at least some of the excited states can be obtained 
by insertion of current operators. 

It is crucial for the construction used above that the 
system is one-dimensional with periodic boundary con¬ 
ditions. This can be seen, for example, noting that the 
ansatz is a decomposition into momentum space 

eigenstates. However, the states obtained by insertion of 
current operators might still describe low energy eigen¬ 
states of an interesting, two-dimensional system with a 
different Hamiltonian. We plan to investigate the prop¬ 
erties of these states in two dimensions in future work. 
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Appendix A: Commutator Between H and i 


In this section we derive an expression for the com¬ 
mutator This commutator was used in Sec¬ 

tion 


HI B|to determine the action of "H on a state ip. 


The explicit form of the Hamiltonian in terms of spin op¬ 
erators iPl 


H — g ^ ^ ( WijWiktjtj, 

k,j 

+ (Ai) 
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We use the Fourier transforms 

Wij Jo, if A; = 0 

^^^j^\{2k-N)h+iN, iffc = l,2,...,iV-l 

r {N - - 2/3ySiN, 

V -v ^ij J if /c = 0 

(^V6 - 2(fc - Nl2f - 2/3)4+iiV, 

[ itt = l,2.iV-1. 

(A2) 


to rewrite "H in terms of . Here (5m = 1 if w mod N = 0 
and 6m = 0 otherwise. A derivation of these formulas is 
given in Section of the Appendix, see also Ref. [551 
where these Fourier sums were evaluated using contour 
integrals. 

We find 


N-l 


\ 1 + 2N‘^ — 9Nk + j, ^ 

k=0 


(A3) 


Starting from this expansion, we next compute the com¬ 
mutator between H and j. Using 

[^mi'^n] — ‘^^abcUn+mi (-^4) 

which directly follows from the commutator algebra of 
the spin operators tf, we obtain 




— ‘2jiSahc^—\^ j^^abc ^ ^ rUZi_m'^m —1 
m—1 

+ i3-N)uy. (A5) 


It is possible to rewrite the sum over Fourier transformed 
spin operators in terms of the operator Cf, which was de¬ 
fined in Eq. (34). To this end, we computed the Fourier 
expansion of C“. The result is 


o ^ Qa 

b E w = + (2p + 2 - 2N)u% 

3 = 1 

2i 

- -^eabc E rnutmuy-p, (A6) 


invariance, 


i{^N) * 


... 

i^j * 3 


E2r = ^^-^E ¥• 


(Bl) 


(B2) 


In order to evaluate the remaining sums, it is useful to 
compute the Fourier sums 


N-l 

k^O 

for n = 1 and n = 2. For j = 0, we have 

N-l 


E^ = 


N-l 


E^^ = 


fc=0 


N{N-l) 


N{N - 1){2N - 1) 
6 ' 


(B3) 


(B4) 


(B5) 


For J = 1, 2,..., A^ — 1, we use the generating function 

-^-1 1 iijN 


f{^) = E 


^iujk _ 


I - 


N-l 


E ^”4 = 


fc=0 


/c=0 

= (- 
V iduj 


I - 6*“^ 

/(w) 


(B6) 


uj—2-Kj/N 


Taking the hrst and the second derivative, we find 
N. 


E = y “ i) 


N-l 


N 


E = y (I - ^ + - w]n). 


(B7) 




Taking the inverse Fourier transforms of these equations 
and solving for W 3 n/z’j and '^jN/zy we 

arrive at Eq. (A2). 


Appendix C: Decoupling Equation for States with 
Current Operators 


for p S {I,..., A^ — I}. Note that for p = I we find the 
same sum over Fourier transformed spin operators that 
occurs in the above expansion of the Hamiltonian (Eq. 
(A5)). Inserting the p = I Fourier mode of Cf , we arrive 


at 


=(N-l)u‘^_,+ieabcuyT^ + lY.-- (A7) 


2 = l 


Appendix B: Fourier Transform of it;-Functions 

The Fourier transforms of Wij = (zi -I- Zj)l{zi — Zj) 
and ref - can both be reduced to one sum by translational 


The starting point for the derivation of the action of 
Cf on a state 

V^afc...ai ('^l; • ■ • 7 ^n') 

= ((/.,, (zi). (z^)(jyi... y\)(o)) (Cl) 

is the null operatoil^ 

{KI^)i{j’Liipsi){zi), with (C2) 

{Kb)^ = ‘^{6ab - isabety. 

When inserting the null operator into a correlation func¬ 
tion of primary fields, the resulting correlator vanishes. 
We therefore have 
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0 = (if,(zi)... ...ci,,AzN){r\rv ■ • ■ ^-\)(0)) 

=i 2^^ ■ ■ ■ i 

dwk 


dwi 


^ ^ 2Triwk Jo 2Triwi J^. 27ri(z — Zi) 


dz 


{4>s-, (zi) . . . J^(z)(j)si (Zi) . . . 4 >sm {zn)J°'’’ {Wk) ■■■ (wi)) 


k 

- wxE j 

9=1“' 


1 

dwk 

r dici 


fo 

2Kiwk 

Jo ‘^KlWl , 

E 27r*(z - 

I 

dwk 

I J 

r dz 


2'Kiwk Jo 2TTiwi 27ri(z — Zi) 


{Js, {zi)... J\z)Js^ {zi) ...J,, (z^) J“'' (wk)... {w ^)). (C3) 


Inserting the OPE between a primary field and a current operator (cf. Eq. @) and the OPE between two current 
operatorJ^, 


r(z)r{w) 


Oab 


2(z — wY 


d~ i^abc 


JYw) 


z — w 


(C4) 


and writing $s(z) = Jsiizi)... 4 >s^{zn), we get 

o = W).E 

j-f Jo 2mwk Jo 2Tnwi J 

k 


dz 


Y 

0 


27ri(z — Zi) z — Zj 


{<^YY)rYwk)--.rjwY) 


9=1 


dwk 


dwi 


2'Kiwk Jo 2.Kiwi 

dz Sba„ 


2Ki{z — Zi) 2(z — Wq)'^ 

k 


($s(z)J“'‘('u;fc) ... .. . J°’Ywi)) 


b 


drcfe 


dwi 


dz i^ban 


y r 

^ Jo 2KiWk Jo 2kiwi J^^ 2ki{z - z^) z - 
dwk 


-{^Y^)rYwk)-.-JJwq)...rY'wi)) 


Y 

0 


Zi - Zj Jo 2Kiwk 


dwk 


dwi 


1 


2Kiwk Jo 2.Kiwi 2{wq — Zi) 


= -(ifn. E 

3{¥^i) 

+ (k;.).E 

q^l 
k 

-ra.E 

9=1 "'0 

fb k /j^a \ 

= E +E 


Jo 2kiwi 

($s(z)J“''(uife) ... J“«+H'u;q+i)J“’-i(w,_i)... J'^Jwi)) 


dwfc 


dwi ieba„ 


2KiWk Jo 2,KiWl Wq — Zi 


{M^)rYwk)---JJwq)...j^Ywi)) 


9=1 


2z2 




-yiKb)J^ba, 

9=1 


diet 


, . r 0— — —{M^)j‘^Ywk)--.jjwq){rrj..r\)m 

ZmWk Jo Wq — Zi 

k q—1 


= -(Kbhy-^Ja,...a,+y „ 2 Jak-’-aq+iaq-i.-.a-i -EE {K)^)Jebag 

9=1 9=ln=-l 


dwk 


din w ^ ^ 

^ ^ (d>,(z)... r^+jwq+YiJ^j-y ■ ■ ■ j-\)m 


2KiWk Jo 2.KiWq Wq — Zi 

fb k /j^a \ 

= E + E 


■f_L-\ ^3 1 


2z2 
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k q—1 

+ E E ^ . . . Jlr J'nJ-V ■ ■ ■ •/-\)(0)) 

q—1 n—— 1 

Multiplying this equation by 2zi gives 


(C5) 




1(5^0 


9=1 


k q—1 


2W)*E E 






($,(z)(j“^i... ... j“\)(o)) 


n——l 


The application of the spin operator T“ to ipak-.-ai can 
be written in terms of the Levi-Civita symbol, 

TVa....a, = {<ilsA^l)...^sA^N){jUl\-.-r\m) 

N 

~ * ^ ' £6a,cV'afc...ag4.ica,_i...(ji • (^6) 

9=1 

Using this result, noting that {K^)it\ = 0 and that the 
operator Cf can be written as C“ = 
we get 




Ci Ipak-.-ai 

A , 

9=1 

+ 2 ra.EE^ 

q^2 n^O 

X ($,(z)(j“^i... ... j“\)(o)). 


is the translation operator and Vij = 2t“t“ + 1/2 is the 
operator that permutes the states of spin i and spin j. 
Note that the momentum is defined modulo 27r. 

The wave function ijjQ as given in Eq. (101 is equivalent 
to 


'0o('5i,...,Siv) (D2) 

N 

= n Y[ (z„ - 

p—1 n<m 

because ipo and ipo differ only by a spin independent con¬ 
stant. 

Applying T to ipo gives 

N 

TAo = <5s n (J33) 

p^l 

N-1 N 

X n n 

n—l m—n-\-l 


Appendix D: Lattice Momentum of Wave Functions 

Following Ref. [SI we define the lattice momentnm op¬ 
erator as p = —iln(T), where 

T = Vn,n-i'Pn-i,n-2 ■ ■ ■'P 2 P (Dl) 


Therefore the momentum of ipo is 


Po 


TT for N/2 odd 
0 for N/2 even 


(D4) 


With Tm“iT ^ = e it follows that •tpak...ai 

has the momentum po ~ 2Trk/N. 


^ A. M. Lauchli, “Numerical Simulations of Frustrated Sys¬ 
tems,” in Introduction to Frustrated Magnetism, Springer 
Series in Solid-State Sciences, Vol. 164, edited by 
C. Lacroix, P. Mendels, and F. Mila (Springer Berlin Hei¬ 
delberg, 2011) pp. 481-511. 

^ N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. 
Teller, and E. Teller, “Equation of State Calculations 
by Fast Computing Machines,” The Journal of Chemical 
Physics 21 , 1087-1092 (1953) 

^ D. C. Handscomb, “The Monte Carlo method in quantum 


statistical mechanics,” Mathematical Proceedings of the 
Cambridge Philosophical Society 58, 594-598 (1962) 

S. R. White, “Density matrix formulation for quantum 
renormalization groups,” Phys. Rev. Lett. 69, 2863 (1992) 

® G. Moore and N. Read, “Nonabelions in the fractional 
quantum hall effect,” Nuclear Physics B 360 , 362 - 396 
(1991) 

° J. 1. Cirac and G. Sierra, “Infinite matrix product states, 
conformal field theory, and the Haldane-Shastry model,” 
Phys. Rev. B 81 , 104431 (2010) 











15 


^ A. E. B. Nielsen, G. Sierra, and J. I. Cirac, “Local models 
of fractional quantum Hall states in lattices and physical 
implementation,” Nat. Commun. 4, 2864 (2013) 

® A. E. B. Nielsen, J. I. Cirac, and G. Sierra, “Quantum 
spin Hamiltonians for the SU{2)k WZW model,” Journal 
of Statistical Mechanics: Theory and Experiment 2011 , 
P11014 (2011) 

” H.-H. Tu, “Projected BCS states and spin Hamiltonians 
for the SO{n)i Wess-Zumino-Witten model,” Phys. Rev. 
B 87, 041103 (2013) 

H.-H. Tu, A. E. B. Nielsen, J. I. Cirac, and G. Sierra, 
“Lattice Laughlin states of bosons and fermions at filling 
fractions 1/g,” New Journal of Physics 16 , 033025 (2014) 
H.-H. Tu, A. E. B. Nielsen, and G. Sierra, “Quantum 
spin models for the SU{n)i Wess-Zumino-Witten model,” 
Nuclear Physics B 886, 328 - 363 (2014) 

R. Bondesan and T. Quella, “Infinite matrix product states 
for long-range SU{N) spin models,” Nuclear Physics B 
886, 483 - 523 (2014) 

E. Ardonne, P. Fendley, and E. Fradkin, “Topological 
order and conformal quantum critical points,” Annals of 
Physics 310 , 493-551 (2004) 

A. E. B. Nielsen, J. 1. Cirac, and G. Sierra, “Laughlin 
Spin-Liquid States on Lattices Obtained from Conformal 
Field Theory,” Phys. Rev. Lett. 108 , 257206 (2012) 

F. D. M. Haldane, “Exact Jastrow-Gutzwiller resonating- 
valence-bond ground state of the spin-1/2 antiferromag¬ 
netic Heisenberg chain with 1/r^ exchange,” Phys. Rev. 
Lett. 60 , 635-638 (1988) 

B. S. Shastry, “Exact solution of an 5 = 1/2 Heisenberg 
antiferromagnetic chain with long-ranged interactions,” 
Phys. Rev. Lett. 60 , 639-642 (1988) 

F. Calogero, “Ground State of a One-Dimensional A^-Body 
System,” Journal of Mathematical Physics 10 , 2197-2200 


F. Calogero, “Solution of the One-Dimensional A^-Body 
Problems with Quadratic and/or Inversely Quadratic Pair 
Potentials,” Journal of Mathematical Physics 12 , 419-436 
(1971) 

B. Sutherland, “Quantum Many-Body Problem in One Di¬ 
mension: Ground State,” Journal of Mathematical Physics 
12 , 246-250 (1971) 


B. Sutherland, “Exact Results for a Quantum Many-Body 
Problem in One Dimension,” Phys. Rev. A 4, 2019-2021 
(1971) 

F. D. M. Haldane, ““Spinon gas” description of the S = 
1/2 Heisenberg chain with inverse-square exchange: Ex¬ 
act spectrum and thermodynamics,” Phys. Rev. Lett. 66, 
1529-1532 (1991) 

““ F. D. M. Haldane, Z. N. C. Ha, J. C. Talstra, D. Bernard, 
and V. Pasquier, “Yangian symmetry of integrable quan¬ 
tum chains with long-range interactions and a new descrip¬ 
tion of states in conformal field theory,” Phys. Rev. Lett. 
69 , 2021-2025 (1992) 

J. C. Talstra, Integrability and Applications of the 
Exactly-Solvable Haldane-Shastry One-Dimensional Quan¬ 
tum Spin Chain, Ph.D. thesis, Princeton University (1995), 
arXiv:cond-mat/9509178. 

P. Bouwknegt, A. W. W. Ludwig, and K. Schoutens, 
“Spinon bases, Yangian symmetry and fermionic represen¬ 
tations of Virasoro characters in conformal field theory,” 
Physics Letters B 338 , 448 - 456 (1994) 

V. G. Knizhnik and A. B. Zamolodchikov, “Current alge¬ 
bra and Wess-Zumino model in two dimensions,” Nuclear 
Physics B 247 , 83-103 (1984) 

P. Di Francesco, P. Mathieu, and D. Senechal, Conformal 
Field Theory, Graduate Texts in Contemporary Physics 
(Springer, 1997). 

D. Gepner and E. Witten, “String theory on group mani¬ 
folds,” Nuclear Physics B 278 , 493 - 549 (1986) 

B. A. Bernevig, D. Giuliano, and R. B. Laughlin, “Coor¬ 
dinate representation of the two-spinon wave function and 
spinon interaction in the Haldane-Shastry model,” Phys. 
Rev. B 64 , 024425 (2001) 

B. S. Shastry, “Taking the square root of the discrete 1/r^ 
model,” Phys. Rev. Lett. 69 , 164-167 (1992) 

K. Schoutens, “Yangian symmetry in conformal field the¬ 
ory,” Physics Letters B 331 , 335 - 341 (1994) 

J. D. Hunter, “Matplotlib: A 2D Graphics Environment,” 
Computing in Science & Engineering 9 , 90-95 (2007) 

P. Ramachandran and G. Varoquaux, “Mayavi: 3D Vi¬ 
sualization of Scientific Data,” Computing in Science & 
Engineering 13 , 40-51 (2011) 









